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Objects driven through periodically modulated potential-energy landscapes in two dimensions can 
become locked in to symmetry-selected directions that are independent of the driving force's orien- 
tation. We investigate this problem in the overdamped limit, and demonstrate that the crossover 
from free-flowing to locked-in transport can depend exponentially on an object's size, with this 
exceptional selectivity emerging from the periodicity of the environment. 



I. INTRODUCTION 

The theme of transport through modulated potential- 
energy landscapes pervades solid-state physics and arises 
in many natural and industrial processes. This problem 
has been studied extensively in the quantum-mechanical 
limit. Considerably less attention has been paid to the 
classical limit, where effects such as viscous damping and 
thermal randomization complicate the analysis. This 
article focuses on noninertial transport of classical ob- 
jects driven through periodically modulated potential- 
energy landscapes by constant, uniform forces. The one- 
dimensional variant of this problem has been thoroughly 
investigated 0, and its results have been applied prof- 
itably to such processes as gel electrophoresis. We focus 
instead on the overdamped motions of classical objects as 
they flow through two-dimensional periodic landscapes, 
about which far less is known. Such higher-dimensional 
periodic landscapes have shown exceptional promise in a 
new category of sorting techniques. Our discussion draws 
upon recent experimental realizations of this process in 
which macromolecules or mesoscopic colloidal particles 
are observed while moving through arrays of microfab- 
ricatedposts and through regular arrays of optical 
traps UELI3- ^ n both cases, particles' differing inter- 
actions with the physical landscape and their differing 
responses to the external driving force can cause them to 
follow radically different paths, thereby providing a novel 
basis for dispersing small fluid-borne objects into distinct 
fractions. 

Section [H] introduces the theoretical framework for de- 
scribing driven objects' interactions with inhomogeneous 
environments in the context of recent experimental real- 
izations. We then apply this in Sec. lIIII to the particularly 
simple case of transport across a linear barrier or poten- 
tial trench. Such a landscape can continuously sort mix- 
tures of objects into two distinct fractions, but with only 
algebraic sensitivity to properties such as size. General- 
izing to periodic landscapes in Sees . II VI and Ivl leads gen- 
erally to fractionation with exponential size selectivity. 
Exploiting this exceptional resolution for practical sepa- 



rations may be difficult, however, in the most straightfor- 
ward implementations. Other potential landscapes, such 
as a line of discrete potential wells, discussed in Sec. IVII 
offer exponential size selectivity with good prospects for 
practical implementations. 



II. MOTIONS THROUGH LANDSCAPES 

A. The Equation of Motion 

Consider a Brownian particle moving, under the influ- 
ence of a uniform driving force Fo, through the force field 
F(r) due to an inhomogeneous medium or landscape. Its 
trajectory is described by the Langevin equation [J, ||| 

^=F(r)+F +r(t), (1) 

where £ is the particle's viscous drag coefficient, and T 
describes random thermal fluctuations. This Langevin 
force satisfies (T(t)) = and (T(t)-T(t+T)) =£k B T6(T) 
at temperature T, where S(t) is the Dirac delta function. 
A sphere of radius a immersed in an unbounded fluid of 
viscosity rj, for example, has £ = Q-Kna. 

In the limit that Fo and F both greatly exceed the scale 
of thermal forces, T, the Langevin equation, Eq. (JTJ, re- 
duces to a first-order deterministic equation of motion. 
This article focuses on two-dimensional systems, the sim- 
plest case exhibiting novel behavior. Even this decep- 
tively simple system yields surprising results, as we will 
see. 



B. The Driving Force 

In the particular case of fluid-borne colloidal particles, 
a uniform driving force might be exerted by viscous drag, 
by gravity, or through electrophoresis, magnetophoresis, 
or thermophoresis. Each of these plays a central role 
in practical fractionation techniques 0. More generally, 
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analogous results should be expected for such related sys- 
tems as electrons flowing through a periodically-gated 
low-mobility two-dimensional electron gas 0, magnetic 
flux quanta creeping through patterned type-II supercon- 
ductors |9L llOt or Josephson junction arrays and 
atoms migrating across crystal surfaces ^jj . 

In some instances of practical interest, the driving force 
itself can be modulated by the physical landscape, lead- 
ing to additional interesting effects 01- These, however, 
are beyond the scope of the present discussion. Time- 
dependent driving forces also lead to exciting new phe- 
nomena, but are not required for the effects we describe. 
We consider the simplest case, where the driving force Fo 
is both uniform and constant and is oriented at a fixed 
angle 9 with respect to the landscape symmetry axis, here 
denoted x. 

In the absence of other influences, particles would 
travel along the driving direction while dispersing diffu- 
sively in the transverse direction. Differential dispersion 
by transverse diffusion has proved useful for continuously 
fractionating heterogeneous samples across laminar flows 
in microfluidic channels |15|. Adding a modulated sub- 
strate opens up new modes of separating particles accord- 
ing to their sizes, and can greatly improve the resolution 
of such separations. 



C. Creating Landscapes 

Several approaches have been introduced in recent 
years for structuring potential energy landscapes on 
molecular, macromolecular, and cellular levels. Among 
these are arrays of lithographically defined microscopic 
posts integrated into hermetically sealed fluidic channels, 
which provide a periodic and precisely tuned alternative 
to the gels used for electrophoresis Arrays of inter- 
digitated electrodes [T^.ll8j also have been used to estab- 
lish periodic potentials through dielectrophoresis. The 
emphasis in these studies, however, has been on ratchet- 
like behavior induced by time-dependent potentials. 

More recently, techniques have been developed for tai- 
loring extensive potential energy landscapes using forces 
exerted by light. The most capable of these exploit op- 
tical gradient forces, meaning that dipole moments in- 
duced in illuminated objects respond to gradients in the 
illumination's electric field. Such forces are the basis for 
the single-beam optical trap known as an optical tweezer 
[l9j . which acts as a potential energy well for particles 
with appropriate optical properties. More generally, an 
extended optical intensity distribution will produce an 
associated potential energy landscape. 

The most straightforward way to project periodic in- 
tensity profiles is to create a standing-wave interference 
pattern from two or more coherent beams of light. Such 
patterns have come to be known as optical lattices, par- 
ticularly when applied to controlling the distributions 
and motions of matter. More general intensity patterns 
can be created with holographic optical tweezers (HOT) 



pfl I2H |22T| or with the generalized phase contrast (GPC) 
technique |23j ] , which establish extended optical trapping 
patterns using computer-generated holograms. 

D. Form Factors 

The physical landscape may be represented by a func- 
tion I(r) describing a potential-determining property 
such as the local optical intensity. An object's potential 
energy at r is determined, not only by /(r), but also by 
the object's response to it. For example, larger particles 
approaching a well-localized optical trap encounter the 
trap's intensity gradients at larger ranges than smaller 
particles. The observation that different objects pass- 
ing through the same environment experience different 
potential-energy landscapes provides the foundation for 
the results that follow. 

The effective potential may be expressed as the convo- 
lution 

F(r) = (/o/)(r) (2) 
= J f(x-r)I(x)d 2 x, (3) 

of the two-dimensional landscape, I(r), with a form fac- 
tor /(r) describing the object's interaction with the land- 
scape. In comparing to experimental realizations, we as- 
sume that contributions from the form factor's third di- 
mension have been integrated out. If J(r) has a symme- 
try axis along the x direction, then the associated force, 

F(r)=-V(/o7)(r), (4) 

generally does as well. Convolving with /(r) broadens 
features in J(r) by an amount that depends on the ob- 
ject's size, shape, orientation, and composition. 

In many cases of practical interest, the convolution in 
Eq. @ is most easily performed using the Fourier con- 
volution theorem: 

(/o/)(r)=^- 1 {/(k)/(k)} (5) 

where /(k) and I(k) are the Fourier transforms of /(r) 
and I(r), respectively, and g(r) = ^ r_1 {g(k)} denotes 
the inverse Fourier transform of g(k). In some particu- 
larly simple cases, both /(k) and 7(k) can be factored 
into components along the x and y directions, reducing 
Eq. (J5J to a product of one-dimensional integrals. In 
other cases, separable approximations for the form fac- 
tor emerge as the leading-order cumulant expansion of 
/»■ 

For example, the form factor for a uniform dielectric 
cube of side a aligned with the x axis and illuminated by 
collimated light is 

/(r)=aae(|-|a:|)e(|-|v|) (6) 
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where Q(x) is the Heaviside step function, and 



a 



2tt 



2e 



(7) 



describes the matter-light interaction, in the quasi-static 
limit, for a material of dielectric constant e immersed in a 
medium of dielectric constant eo [24[ . Both this geometry 
and the collimated light field are far simpler than would 
be encountered in most real-world optical trapping imple- 
mentations pl|. but serve to illustrate our approach. A 
more complete treatment of optical forces also would in- 
corporate polarization effects, which cannot be captured 
in the present scalar theory. Higher-order effects such as 
Mie resonances |U could be taken into account through 
a, but will be ignored in the current discussion. Note that 
a is negative for a high-dielectric-constant material in a 
low-dielectric-constant medium; such particles are drawn 
toward regions of high intensity. Low-dielectric-constant 
particles, by contrast, are repelled by light. 

The aligned cube's form factor is separable, with 
Fourier transform 

/(k) = aa 3 f x (k x a) f y (k y a). (8) 

The individual components are readily shown to be 

sin ka 



fx(ka) = fy(ka) = 



ka 



Their leading-order cumulant expansion 

1 



fx{ka) = f y {ka) w exp ( -- k a' 



(9) 



(10) 



for ka < 7r, demonstrates that the form factor's Fourier 
transform depends sensitively on particle size for a given 
wavenumber. Note that, as defined, f x (ka) and f y (ka) 
are dimensionless and normalized to unity at ka = 0. 

The form factor for a uniform dielectric sphere of 
radius a illuminated by collimated light of wavelength 
A > a is HI, 



/(r) = a V a 2 — r 2 0(a — r), 



(11) 



which is not separable. The leading-order cumulant ex- 
pansion of /(k), however, is separable, with 

f(k)^a^-f x (k x a)f y (k y a), (12) 

where 



Mka) = fy(ka) = exp f -— fcVJ , (13) 
for ka < 8. 

More generally, an object's form factor is nonzero only 
over a limited domain, set by its size. The corresponding 
Fourier transform thus depends strongly on ka, within 
the appropriate range of wavenumbers. We capture the 



ramifications of this boundedness by adopting the sepa- 
rable Gaussian form 



(14) 



whose Fourier transform 

/(k) = 27raa 3 f x (k x a) f y {k y a) (15) 
has components 

1 



f x (ka) = f v {ka) = exp ( -^-k 2 a 2 



(16) 



Which wavenumbers come into play depends on the land- 
scape, I(v). The following Sections explore a few partic- 
ularly effective choices. 



III. LINEAR FRINGES 

In part to motivate a discussion of periodic potential 
energy landscapes, we first consider how objects traverse 
a single trench or barrier arranged at an angle to the 
driving force. This kind of landscape may be realized, 
for example, by creating a linear optical trap with a 
cylindrical lens or a diffractivc line generator. Because 
such an optical landscape can act as either a barrier or a 
trench, depending on the sign of a, we will refer to both 
as fringes. In either case, a fringe aligned with the x axis 
inhibits transport in the transverse direction. We model 
the landscape as a Gaussian profile of intrinsic width w, 



I(r) = I exp 



y 

2w 2 



(17) 



Using the object's form factor as defined in Eq. (|14|) . the 
associated potential is 



V(r) = 2nal 



cr(a) 



exp 



2a 2 (a) 



(18) 



The fringe's apparent width to a particle of size a is 
broadened to cr(a) = \/ a 2 + w 2 . 

In the limit that thermal forces may be ignored, the 
equations of motion reduce to the deterministic form 



dx 
II 
dy 
dt 



v cost) 

r 1 F y (y)+v sin I 



(19) 
(20) 



where the landscape- free drift speed is vq = £ 1 Fo and 

,,2 



Fy{y)=2nal { 



o 



■ V exp 



y 



2a 2 (a) 



(21) 



The landscape's restoring force F y (y) reaches a maxi- 
mum at a distance y = y max from the fringe's axis, with 
|/max = c(a) for our particular example. If Fosin# > 
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FyiVmax) then a particle can cross the barrier. Such par- 
ticles may be said to escape the barrier. By contrast, 
particles for whom the barrier is insurmountable travel 
unimpeded along the x direction at speed v x = vq cos 9. 
Such particles are said to be locked in to the landscape. 

The marginal angle 9 m at which an object just barely 
remains locked in to the barrier determines which ob- 
jects are deflected and which are not. The dependence 
of 6 m on particle size and other characteristics estab- 
lishes the sensitivity of the sorting technique. Referring 
to Eqs. H2U|) and (|21|) . the condition for locked- in trans- 
port, 



sm < sm m = 

\a\I 2-7T a 3 w 



Fq ^/e a 2 + w 2 



(22) 
(23) 



applies both to attractive trenches (jj max = +er) and 
repulsive barriers (y max = — c)- The general result, 
Eq. i|22|) . applies even if /(r) is not separable because, 
in this case at least, I(r) is independent of x. 

The particular result in Eq. (|23|1 shows that the 
marginal lock-in angle depends only algebraically on size, 
and only linearly on other properties through a. This is 
neither better nor worse that the performance offered by 
other established techniques such as gel electrophoresis 
or flow- field fractionation 7] . One substantial benefit of- 
fered by selective transport across a fringe is its ability 
to process a continuous stream of objects rather than be- 
ing restricted to discrete batches. The selected fraction, 
moreover, can be tuned continuously, for example by ad- 
justing Jo, Fa, w, or 9. Optical implementations also can 
be optimized by varying the wavelength of light, in which 
case resonances might be exploited as a complementary 
mechanism for size separation. 

Although a single fringe's performance is somewhat 
lackluster, one might expect multiple fringes to fare bet- 
ter. The first step along this direction is to consider a 
pair of parallel Gaussian fringes. The effective potential 
is the sum of two single-fringe potentials: 



V(r) = 2-itaI 



a (a) 



exp 



2a 2 



exp 



(y-b/2Y 

2a 2 



(24) 



where b is the fringe separation. If b < <r, the two fringes 
overlap enough that the landscape resembles a single, 
broadened fringe. Again, no more than algebraic selec- 
tivity should be expected. In the opposite limit, b 3> er, 
the fringes are independent, and particles cross the dou- 
ble barrier with the same facility with which they cross 
one. Neither of these cases offers benefits over the single 
fringe. 

For intermediate fe, on the other hand, the landscape 
consists of two unequal barriers, the smaller of which lies 
between the two fringes. The smaller barrier's height 



depends strongly on b/a, which, in turn, depends on 
the particle size a. This lower intermediate barrier does 
not affect the fringes' overall ability to separate objects, 
which is dominated by the larger barrier. It suggests 
though possibility, that transport across N overlapping 
fringes could be highly sensitive to particle size. Those 
particles not able to jump the inter-fringe barriers will be 
locked in and swept aside while others will hop from one 
fringe to the next across the field. Highly selective sort- 
ing is thus possible if edge effects due to the first or last 
fringe (in the case of trenches or barriers, respectively) 
can be circumvented. 



IV. SINUSOIDAL LANDSCAPES 




FIG. 1: Schematic representation of a sinusoidal landscape, 
modulated along the y direction. 

The foregoing discussion suggests that a periodically 
modulated landscape inclined at an angle to the driving 
force might be more effective than a single fringe at sort- 
ing objects by size. To make this more concrete, and to 
illuminate the role of periodicity, we consider the sim- 
plest and most instructive example of such a landscape, 
a sinusoid in the y direction: 



J(r) = I Q cos(fc y), 



(25) 



shown schematically in Fig. ^ Apart from its mathemat- 
ical simplicity, this landscape has the advantage of being 
readily implemented experimentally. In the case of opti- 
cal forces, a sinusoidal pattern can be created by interfer- 
ing two coherent laser beams, with the spatial wavenum- 
ber fco determined by the optical wavelength and the an- 
gle between the beams. Such an interference pattern 
is known as a one-dimensional optical lattice |26(, and 
is commonly used to control and distribute cold atoms. 
More recently, optical lattices have been used to separate 
populations of colloidal particles on the basis of their sizes 
and indices of refraction Q . 

The key to such a potential's efficacy is in its Fourier 
transform: 



7(k) = (2ir) 2 I 6{k x )5(k y -k ). 



(26) 
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Convolution according to equation [5] then results in the 
potential landscape picking out the component of the 
form factor's Fourier transform at wavenumber k^\ 

V(r) = I f(k)S(k x )S(k y - k ) cos(k x x) cos(k y y) d 2 k 



(27) 

= Jo/(0,Ao) cos(fcoy). (28) 
Assuming a separable form for /(k) such as Eq. JTBJ, 
V(r) =I f v (k a) cos(fe y) (29) 
= a\/27r/o a 3 exp y~ & 2fl2 ^ cos(fc y). (30) 

For particle sizes such that kn_a > a/3j the amplitude 
of the landscape's sinusoidal modulations now depends 
strongly on particle size through the wavenumber depen- 
dence of the form factor. The particular form in Eq. i|3U|) 
reflects our choice of a Gaussian form factor in Eq. (|14|l . 
However, the arguments leading to this choice reveal that 
comparable results should be obtained quite generally for 
particles whose size is smaller than the wavelength of the 
physical landscape's undulations (e.g., for spherical par- 
ticles, for koa < 8). 

A. Deterministic Limit 

A particle's trajectory through the sinusoidal land- 
scape is described in the deterministic limit by Eqs. i|19|) 
and J3UJ, with 



F y (y) = av27r/ofco« exp 



1 



7 2 2 



sin(fc y). (31) 



Non-separable form factors again yield similar results, 
because I(r) is independent of x. Because sm(k y) < 1, 
particles become locked into the x direction for orienta- 
tions satisfying 



sin 9 < sin 9 r , 



\a\Io 

Fn 



2-Kkoa exp 



kna 2 



(32) 



By contrast to Eq. 1)23(1 . this reflects an exceptional sort- 
ing sensitivity: whether or not a particle becomes en- 
trained by the fringes depends exponentially on parti- 
cle size for kna > 1. Among established fractionation 
schemes, only affinity chromatography offers comparable 
selectivity Q , and this can operate only on discrete sam- 
ples of a limited class of macromolecules. Fractionation 
in a sinusoidal landscape, by contrast, can operate on 
continuous sample streams and can be implemented for 
a wide range of sample types. 

Equations i|T$|) . (|2T)jl . and (I5U|) can be directly inte- 
grated for this simple landscape. The motion in the x 
direction is trivial: 



x(t) = v cos9t. 



(33) 



If the particles are locked in (i.e., if Eq. I|32|) is satisfied), 
then the particles make no progress in the y direction, 
and y(t) is constant in steady state. Otherwise, integra- 
tion gives 



y(t) = — arctan 
k 



i sin ( 



V 



sin t 



V 



( k O y ot l. 2a 2 

tan I — - — y sin — n z 



(34) 

where the relative strength of the landscape's modulation 
is 



77(a) 



I k o f(0, k a) 

Fa 



(35) 



This motion can be seen by inspection to be periodic, 
with the period 



T = 



2tt 



knVo\/ sin 2 6 — rf 



(36) 



corresponding to the time to advance one fringe spacing, 
b = 27r/fco. The mean velocity in the y direction is thus 
given by 



v \ sin 



(37) 



On average, the particle travels at an angle to the x 
axis, given by 

(v„) f 0, sin6»<r/ 

tan ^H = WiZZEZ (38) 
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FIG. 2: (a) Travel direction as a function of orientation for 
an inclined sinusoidal landscape as a function of orientation 
for fixed size and 77(a) = 0.1, 0.2, and 0.3. Trajectories are 
locked in to tp(9) — for 6 < 8 m . The diagonal dashed line 
indicates the result with no landscape, (b) Deflection angle 
as a function of particle size a at fixed driving orientation 
tan# = 0.441, assuming 770 = 0.4, independent of a. 

The deflection angle is plotted in Fig. OJa) as a func- 
tion of the of the driving force's orientation 9 for various 
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values of the normalized potential 77(a). The direction in 
which particles flow increases from tp = as the driv- 
ing force's orientation crosses the condition for marginal 
lock-in, 9 m = arcsinf7. At steeper angles, tf> approaches 
9. 

While this result is quite general, we can make the 
dependence on particle size more explicit by assuming the 
following functional form for 77(a), implied by Eqs. i|14|) 
and (35} : 



77(a) = 770 exp ( -\k t 



(39) 



where 770 = 2-Kaa 3 Ioko/F . Fig.^b) shows the resulting 
dependence of deflection angle on particle size, if we as- 
sume that the driving and trapping forces are adjusted 
such that 770 is a constant, independent of a. It can be 
seen that particles which are not locked in to the fringes 
at ip(a) = are fanned out into various directions, de- 
pending on their size. 

Unlike the case of the single fringe, where a particle 
either flows along the fringe or else travels in the driv- 
ing direction, the sinusoidal landscape's continuous dis- 
persion distributes heterogeneous samples into multiple 
fractions, but also limits the achievable size resolution. 
The fraction dispersed into a finite angular range At/j 
around ?/> includes an associated range of sizes 



Aa w I — 

oa 



AV>, 



(40) 



which, for the locked-in fraction at i/j = Aip/2 <C 1 is 



A„ = ,-,.s- H { ^ ) At/* 2 
oa 



(41) 



Thus, the exponential size selectivity implied by Eq. i|32|) 
can be lost in the exponentially wide collection window 
imposed by 77(a) on practical implementations. This per- 
formance cannot be improved by passing the set of par- 
ticles through the fringes a second time, because of the 
fixed relationship between Aa and Atp. 

Although single-stage fractionation by a sinusoidal 
landscape yields broad size distributions, a narrow range 
of particle sizes still can be captured by using the fol- 
lowing, two-step process. The deflection angle is first set 
such that all particles larger than a certain size 02 will 
be locked in. These locked-in particles are discarded, and 
the remaining particles are sent through a second poten- 
tial landscape, with a different deflection angle, chosen 
such that all particles larger than a second size ai < a2 
are locked in. Only the locked-in particles from this sec- 
ond stage are then retained, so that all of the remaining 
particles have sizes in the range [01,02], which can be 
made as small as desired. 

With periodicity providing the essential ingredient for 
achieving exponential size selectivity, it might be ex- 
pected that any periodic landscape would do. Unfortu- 
nately, this is not necessarily so. We already have demon- 
strated in Sec. IHII that an array of well-separated Gaus- 
sian fringes offers only algebraic, rather than exponential, 



size selectivity. A more general periodic landscape with 
wavelength 2n/ko can be expanded as a Fourier series: 



(42) 



71=0 



with Fourier coefficients f3 n . If one of these coefficients is 
significantly larger than all the others, then the equations 
of motion can be approximated by Eqs. (T$)) and J2U|, 
and equivalently good size selectivity will be obtained. 
If, on the other hand, no single component dominates, 
then the superposition will not necessarily perform so 
well. 



B. Biased Diffusion 

Modelling thermal effects is reasonably straightforward 
for the sinusoidal landscape. In this case, the Langevin 
equation (Eq. is most readily solved by transforming 
it into a Fokker-Planck equation for the probability den- 
sity p(r, t) of finding particles at position r at time t. If 
inertial effects are negligible, the Fokker-Planck equation 
for motion transverse to the fringes reduces to a Smolu- 
chowski equation pj, 



d tP (y,t) + d y S(y,t) = 0, 



(43) 



where the probability current is 

S(y, t) = r 1 (F(r) - k B Td y ) p(y, t) . (44) 

In this equation, F(r) is the total force on the particle, 
including the driving force and the force due to the po- 
tential landscape, T is the temperature, and ks is Boltz- 
mann's constant. 

Following Ref. 0, Eq. O can be solved in the steady- 
state limit by taking S(y,t) — S(y), independent of t. 
The resulting average drift velocity (v y ) for the sinusoidal 
potential of Eq. (|30|l is given by 



Ml 



Do sin ( 



2sin0 

1 H Im 

V 



Si T, 



sin 9 



V 



(45) 



where, as before, 77 = fco^o/(0, koa)/Fo, and we have 
introduced the normalized temperature t = k^F '/Vq. 
The function Si(t,x) is defined recursively in terms of 
a continued-fraction expansion, 



S„ (t, x) 



1/4 



+ in x + S n+ i (r, x) ' 



(46) 



which converges rapidly with increasing order n. 

The average velocity in the x direction is unchanged 
from the zero-temperature case. The mean deflection 
angle ijj is thus given by 



. , 2 sin 9 
tan ip = tan 9 1 H lm 



Si 



sine 



(47) 
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V. SEPARABLE TWO-DIMENSIONAL 
LANDSCAPES 




FIG. 3: (a) Deflection as a function of orientation at finite 
temperatures, r = 0.01, 0.1, 0.2 and 0.3, assuming r\ = 0.4. 
The diagonal dashed line indicates the result with no land- 
scape, (b) Dependence of the travel direction on particle size 
a for rjo = 0.4, tan9 — 0.441 and the same set of tempera- 
tures. Raising the temperature weakens the size dependence 
of tp(a), and thus reduces the selectivity. 



This is plotted in Fig. GH a ) as a function of the angle 
6 of the driving force, for a fixed value of the normal- 
ized potential 77(a), and for various values of the nor- 
malized temperature r. It can be seen that the effect 
of increasing temperature is to smooth out the transition 
between the locked-in and freely flowing states of motion. 
In the zero-temperature limit, the deflection angle is zero 
for all angles 9 < arcsinf7. For finite temperatures, the 
mean deflection angle is non-zero even in this "locked-in" 
state: the particles have a finite probability per unit time 
of being driven over the inter-fringe barrier by thermal 
fluctuations, and thereby advancing in the y direction. 

The benefits of operating in the deterministic regime, 
in which thermal forces are negligible, can be seen in 
Fig. IHJb) , where the deflection angle is shown as a func- 
tion of the particle size a, for a fixed orientation 0, over 
a range of temperatures. Contrary to previous assertions 
that thermal effects can enhance size selectivity Q, the 
only effect of thermally-assisted hopping in this system 
is to diminish the sorting resolution. 

In other words, achieving high-sensitivity sorting in 
practice will require that the thermal energy scale be 
small compared to the landscape's modulation. In a 
real-world implementation, increasing the depth of mod- 
ulation Iq of the physical landscape will often be more 
practical than decreasing the temperature. Retaining the 
same lock-in conditions then requires that the driving 
force Fq to be increased proportionately. The practical 
limit on the achievable sorting efficiency will then be set 
by the maximum driving force or depth of modulation 
that can be obtained. For example, the limitation for 
an optically-implemented landscape is the available laser 
power. 



FIG. 4: Schematic representation of a landscape sinusoidally 
modulated in both the x and y directions according to 
Eq. QHJ. 

The one-dimensional potential of Eq. i'2b\) is one of 
the few landscapes that allows for exact solutions of the 
equations of motion. Two-dimensional landscapes can be 
solved analytically only if the potential can be written as 
a sum of modulations in the x and y directions. In par- 
ticular, we can consider separate sinusoidal modulation 
with the same period in the two directions: 

Z(r) = I [sin(fc y) + sin(fc j:)] , (48) 

shown schematically in Fig. ^ This landscape is inter- 
esting mainly because it leads to decoupled equations of 
motion: 

^- = C^oh cos(fc x) + v cos 9 + j(t) (49) 
at 

% = rVofco cos(fc 2/) + v sm9 + j(t), (50) 
at 

where ^(t) = r/£. Nevertheless, such a landscape could 
be implemented experimentally using optical forces. For 
example, mutually incoherent pairs of laser beams inter- 
secting at right angles would lead to a potential of the 
form given in Eq. (|48J) . as would pairs with orthogonal 
polarization. 

Since the motions in the x and y directions are inde- 
pendent in this case, the same exponential sensitivity to 
particle size, Eq. IMl'l) is obtained, in the absence of ther- 
mal forces. As well, the same integration can be used 
to determine the average deflection angle for free-flowing 
particles, analogous to Eq. 



w= v^2E2. (5i) 

y cos 2 9 — t] 2 

Similarly, for finite temperatures, the mean deflection an- 
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lie is given by 



tan if) = tan 9 



n 


2 sin £ 


-Im 


5i 






') 








1 4 


2 cos £ 


-Im 


Si 






n 









(52) 




FIG. 5: (a) Deflection as a function of orientation for a sep- 
arable two-dimensionally modulated landscape at 77(a) = 0.1, 
0.2, and 0.3. The diagonal dashed line indicates the result 
with no landscape, (b) Size dependence of the deflection an- 
gle for 770 = 0.4 and tan 6 = 0.441. 

The zero temperature deflection angle is plotted in 
Fig-EJas a function of the angle 9 of the driving force, for 
a fixed value of the normalized potential 77. Also plotted 
is if) as a function of particle size a for a fixed 9. The 
results can be seen to be similar to those obtained with 
one-dimensional fringes. In other words, no qualitative 
difference is obtained in this case by modulating in two 
directions rather than just one. 

In order to see new effects of increased dimensionality, 
it is necessary to consider landscapes that cannot be sep- 
arated into one-dimensional terms; i.e., landscapes where 
the motion in one dimension depends on the position in 
the other. Analytical solutions are not available for such 
landscapes. However, it is possible to develop limiting 
arguments that illustrate the novel features of transport 
in such landscapes, including the continued possibility 
for sorting that is exponentially sensitive to particle size. 



VI. LINEAR TRAP ARRAYS 



A. Periodically modulated fringe 



Combining aspects of Sees. IIIII and IIVI we next con- 
sider landscapes that are uniform outside a bounded re- 
gion in the y direction, and are periodically modulated 
in the x direction. This features the clean separations of 
the continuous barrier with the exponential selectivity of 
sinusoidal landscapes. It also provides a straightforward 
example of the surprising complexity of non-separable 
landscapes. 



The simplest exemplar is the modulated line, 

1 + s cos{kox) 



I(r)=I A(y) 



1 + a 



(53) 



where A{y) describes the transverse profile and is peaked 
at ^4(0) = 1. Here, the factor s controls the depth of 
modulation along the line and falls in the range < s < 
1. Such an array can be realized, for example, with a 
linear array of discrete optical tweezers. Choosing s > 
1 would correspond to alternating potential wells and 
barriers, which also could be implemented optically, for 
instance with two different wavelengths of light. 

At a given driving orientation 9, objects are either 
locked in to the array and deflected, or else escape into 
the driving direction. This is unlike the sinusoidal land- 
scape, for which even the particles that are not locked in 
are deflected away from the driving direction. Collection 
of the desired fraction should thus be more straightfor- 
ward for the linear array of traps. 

The equations of motion for objects driven determin- 
istically through I(r) at angle 9 are 



dx 

~dt 
dy 
dt 



£ 1 F x (r) + v Q cos6> 
C l F y (r)+v sin 9, 



(54) 
(55) 



where the components of the substrate-mediated force 
reflect a convolution with a particle's form factor. We 
will assume the particles' form factors to be separable, as 
in Eq. I|14|l . thereby sacrificing some generality in favor 
of clarity, so that 



F x (r) = aI k a 2 A(y) 



sf x (k Q a) sin(fc :r) 
1 + s 



(56) 



F y (r) = -aha 1 d y A{y) ^— , (57) 

where the effective transverse profile is A(y) = (f y oA)(y). 
Even this simplified set of coupled equations is highly 
nonlinear and cannot be integrated directly. Instead, we 
resort to limiting arguments to determine when particles 
become locked in and when they escape. These estimates 
provide the basis for our claim that sorting by inclined 
arrays of traps or barriers can offer exponential size se- 
lectivity. 

As for the uniform fringe in Sec. lllll the restoring force 
F y (r) attains its maximum value along y — y max , with 
2/max > for attractive wells and y max < for repulsive 
barriers. For this separable model, y max is a solution to 
dyA(y) = 0. In more general non-separable systems, it 
varies with position x along the array. In either case, this 
threshold depends on the object's geometry and compo- 
sition through the form factor /(r). 



A particle's trajectory must cross y 



if it is 



to escape the line of traps. This requires there to be 
at least some points x along the array where Fq sin 9 > 
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F y {x,y laayi ). Limits on this condition can easily be es- 
tablished, with 

F sin 9 > max{F y (x, 7/ max )} (58) 

X 

ensuring that every trajectory escapes and 

F sin6» > mm{F y (x, y max )} (59) 

X 

opening up the possibility that at least some trajectories 
might. The associated bounds on 9 m , the marginally 
locked-in angle, are 

Via) — < sin ro < 77(a) — , (60) 

1 + s 1 + s 

where 77(a) = ctkoa 2 A(y max ) / Fq . With this definition, 
r?(a) depends only weakly on a, as in Sec. Mil 

The limit of weak modulation, s = 0, once again yields 
Eq. (|23|l . the result for a continuous barrier or trench. 
Similarly, since linia^oo / x (fco a ) — 0, large particles with 
k$a > 1 are not significantly affected by the modulation. 
Smaller particles encountering a deeply modulated line, 
s — ► 1, are more interesting. Unfortunately, the simple 
bounds in Eq. (|60|l have no predictive power in this range, 
because lim a ^o fx(koa) = 1 and Eq. (|60[) reduces to < 
sin(9 m < 77(a). 

This is not to say that exponential selectivity is lost 
in this range, but rather that a more detailed analysis is 
required to ascertain when it can be attained. To illus- 
trate the possibility of achieving exponential sensitivity, 
we consider an experimentally realizable landscape con- 
sisting of a line of discrete optical tweezers Q , which we 
model as a line of discrete Gaussian wells. 



B. Line of Gaussian Wells 

After convolution with a Gaussian form factor, a single 
well of intrinsic width w takes the form 

^(r) = -V exp (-£2) (61) 

with a 2 (a) = w 2 +a 2 . This should not be mistaken for an 
accurate model of an optical tweezer's potential well, but 
rather as a tractable model whose behavior approximates 
that observed in actual optical traps. A line of such wells 
separated by a distance 6 results in the potential energy 
landscape 

"W-wl«(-^). m 

Any trajectory locked in to this periodic landscape will 
itself be periodic in x. This means that such a trajectory 
passes through a sequence of turning points at which 
d x y(x,t) = dty(x,t) = 0. Any trajectory lacking such 
turning points cannot be locked in, and so must escape 



from the line of potential wells. Turning points come in 
two varieties: those where particles make their nearest 
approaches to the wells' centers, and those correspond- 
ing to their furthest excursions from the line of traps. 
Particles can escape when the latter type disappear. 

For small to moderate driving angles 9, the more dis- 
tant turning points occur near the midplanes between 
traps, where the restoring force is weakest. Considering 
the influence of just two traps (appropriate for 6 > a), 
centered at x = and x = b, this suggests the point 
of escape will be near x = b/2 and y = a. Expanding 
around this point yields 

sin0 m < 77(a) exp ^-^2) ' ( 63 ) 

where 77(a) = (2/^/e)Vo/(aF ) measures the traps' 
strength relative to the driving force. Because Eq. (|63|l is 
an upper bound, no locked-in trajectories can occur for 
6 > 9 m . Equation therefore establishes the expo- 
nential size dependence of particles' deflection. 

Figure [51 shows results of numerical simulations of 
transport across a line of Gaussian potential wells. These 
simulations were designed to model the experimental de- 
sign of Ref. in which colloidal spheres are driven 
by flowing fluid past an inclined line of discrete optical 
traps. The driving force for this system is Fq ss £au, 
where £ is viscous drag coefficient corrected for hydro- 
dynamic coupling to walls, a is the radius, and u is the 
flow speed. The sample trajectories in Fig. HJa) were 
calculated for w — 0.4 6 and 77(a) = 9.7 a 2 /a 2 at a fixed 
orientation of 6 = 17.5°. The demonstrate that spheres 
with radii larger than a — 0.1 b are locked in to the array 
of traps under these conditions, while smaller spheres 
escape. Even a comparatively short array can resolve 
differences in radius of just a few percent, suggesting 
that nanometer-scale resolution should be attainable for 
hundred-nanometer-scale spheres in practical optical im- 
plementations. 

Figure |5Ib) shows how the marginally locked-in an- 
gle varies with size for this array. The lower dashed 
curve in Fig. H3 is the prediction of Eq. (55} ■ Its very 
good agreement with simulation results in this parame- 
ter range confirms that the limiting argument establishes 
a useful lower bound on the 9 m . These results therefore 
confirm that fractionation by a line of traps offers excep- 
tional size selectivity in an appropriate range of condi- 
tions. Figure E£b) also demonstrates that the locked-in 
fraction can be deflected to large angles, contrary to as- 
sertions in previous reports Q. 

While b > a ensures optical fractionation's exponen- 
tial size selectivity, other considerations provide a ba- 
sis for optimizing the inter-trap separation. The total 
lateral deflection for a captured particle in an A-trap 
array is (N — 1)6 sin#. The array's efficiency can be 
defined accordingly as the lateral deflection per trap: 
A(a, 6) = 6sin6'. Choosing 6 = 2er(a) optimizes this ef- 
ficiency at A = (4/e)Vo/Fj. This result, however, does 
not necessarily optimize sensitivity to particle size. 
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(a) 

a = 0.096/7 
a = 0.098/7 



a = 0.100 fc, 



with 




0.05 0.1 

alb 



0.15 



FIG. 6: (a) Trajectories calculated according to Eqs. I|54|l and 
155t for the line of Gaussian wells described by Eq. 162H . The 
wells are separated by distance b and have intrinsic width 
w = 0.4 6. Their effective width is a — \fw 2 + a 2 , where a is 
the radius of a sphere flowing through the array. The effective 
potential well depth is r/(a) = (2/Ve) Vo/(aF ) = 9.7 a 2 /a 2 . 
With the driving force oriented at 9 = 17.5°, spheres with 
radii larger than a = 0.1 b are locked into the line, (b) De- 
pendence of the marginally locked-in deflection angle 8 m on 
radius, a. The lower dashed curve is the prediction of Eq. 163H 
and the upper from Eq. 1601 . The dashed line indicates the 
orientation along which the data in (a) were calculated. 



The sensitivity may be formulated as 

dA(a,b) 



S(a,b) 



and is optimized by setting 

dS(a, b) _ d 2 A{a,b) 



db 



dbda 



= 0. 



(64) 



(65) 



This yields an optimal separation somewhat larger than 
that for maximum deflection: 



^-2=1 + x(o) + V* + X 2 (a), 



(66) 



X(a) 



if (a) cr(a) 



77(a) cr'(a) 



(67) 



Although fractionation by a line of optical traps has been 
demonstrated in practice Q , optimization based on these 
criteria has yet to be implemented. 

VII. CONCLUSIONS 

Periodic potential-energy landscapes have exceptional 
promise for sorting continuous streams of mesoscopic 
objects. Whether an object becomes locked in to a 
symmetry-selected direction through the landscape or in- 
stead follows the direction of the driving force can depend 
sensitively on size. This can be shown quite generally 
for the separable potentials considered in Sees. IIVI and 
IVl More subtle landscapes, which involve coupled mo- 
tions in two or more dimensions, are more difficult to 
analyze. Approximate arguments and simulations show 
that a particular one of these, a line of Gaussian wells, 
offers both exponential size selectivity and clean binary 
separations. More sophisticated, non-separable, higher- 
dimensional landscapes, such as two-dimensional arrays 
of optical traps 0, optical lattices Q, and microfabri- 
cated post arrays |2| , can distribute continuous distribu- 
tions of objects into discrete fractions j2?|]. The analy- 
sis in this case is made far more difficult by the lack of 
closed-form solutions to the equations of motion, even in 
the deterministic limit. 

Randomization by thermal forces substantially de- 
grades the selectivity with which a one-dimensionally 
modulated landscape can retain objects. A related study 
demonstrates that thermal forcing restructures the pat- 
tern locked-in states in a two-dimensional array of poten- 
tial wells [27T | , and eventually wipes them out as the array 
grows in size. This contradicts the assertion |4( that ther- 
mally assisted hopping can lead to exponential size selec- 
tivity. Fortunately, the sorting processes discussed here, 
as well as their generalizations, can be driven into the de- 
terministic limit by increasing the driving and trapping 
forces. 

Continuous, continuously tuned chromatographic size 
separations should have many applications in biological 
research, drug discovery, and purification of mesoscale 
materials. This Article outlines the basic principles by 
which they work, and suggests considerations for their 
optimization for particular applications. 

This work was supported by the National Science 
Foundation under Grant Number DBI-0233971 and 
Grant Number DMR-0304906. 
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